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NUMERICAL  SIMULATION  OF  THE  LASER-TARGET  INTERACTION 
AND  BLAST  WAVE  FORMATION  IN  THE  DNA/NRL 
LASER  EXPERIMENT 

I.  INTRODUCTION 

Over  the  past  several  years  the  Plasma  Physics  Division  at  the  Naval 
Reasearch  Laboratory  has  developed  an  experimental  apparatus  which  is 
intended  to  simulate  certain  physical  processes  characteristic  of  a  high 
altitude  nuclear  explosion  (HANE)1.  This  program  is  under  the  auspices  of 
the  Division  of  Atmospheric  Effects  of  the  Defense  Nuclear  Agency.  It  has 
the  primary  goal  of  a  physical  understanding  and  reliable  prediction  of  the 
degradation  of  radar  and  communication  systems  through  a  HANE-disturbed 
environment. 

The  first  series  of  experiments  (ending  March  1985)  were  initiated  by 

irradiating  one  side  of  a  thin  foil  target  with  a  neodymium  laser.  Dark- 

field  shadowgraphy  for  the  case  of  a  high  pressure  (  >  1  Torr)  ambient  gas 

taken  at  50  to  150  nsec  after  the  laser  pulse  show  that  the  rapid  heating  and 

evaporation  of  the  target  material  is  followed  by  a  blast  wave  expanding 
2  R 

toward  the  laser  ,  here  denoted  as  the  forward  direction.  Except  for  the 
pronounced  aneurism,  which  is  found  at  high  pressures  (  -5  Torr  )  and 

generally  occurs  near  the  laser  axis,  the  blast  wave  appears  quite  smooth. 
On  the  other  hand  for  those  targets  supported  by  a  stalk,  shadowgraphy 
reveals  that  the  rearward  moving  blast  wave  is  very  disrupted  with  many  small 
flocculi.  The  forward  aneurism  and  rearward  floccuii  are  relevant  tc  the 
question  of  whether  structure  can  be  seeded  in  the  early  stages  of  blast  wave 
Manuscript  approved  January  29.  1986. 


formation.  To  understand  the  mechanisms  responsible  for  these  observed 
features  it  is  clearly  necessary  to  have  a  realistic  model  for  the  initial 
laser- target  interaction  and  the  subsequent  formation  of  the  blast  wave. 
This  report  presents  the  results  of  a  numerical  simulation  of  these  different 
stages.  The  code  includes  a  number  of  pertinent  physical  processes  in  order 
to  provide  a  complete  qualitative  picture  of  the  dynamics.  Existing  measured 
data  on  the  electron  number  density  and  spectroscopy  can  be  used  for 
quantitative  verification  of  the  numerical  simulation.  Since  the  code  is  one 
dimensional  it  cannot  directly  answer  the  question  of  early  time  stucture, 
but  it  is  a  step  toward  a  global  model  for  the  DNA/NRL  laser  experiment  and, 
as  such,  can  be  used  as  a  framework  for  more  detailed  studies. 

The  results  from  the  numerical  models  show  a  hot  (  Tq  -30-70  ev),  dense 
(na  -few  xIO  cm  J  )  cavity  pushing  on  the  forward  moving  blast  wave  after 
-50  nsec.  At  this  time  the  engulfed  ambient  mass  is  much  larger  than  the 
amount  of  ablated  material.  The  gross  structure  of  this  forward  half  of  the 
biowoff  is  similar  to  an  isothermal,  Taylor-Sedov  blast  wave  and  the  outflow 
within  the  cavity  is  subsonic.  This  cavity  however,  i3  filled  with  A £  , 
Ai  ,  A 4  ,  N  ,  and  N  ,  indicating  a  highly  non-equilibrium  plasma.  On 
the  rear  side,  the  blast  wave  i3  colder  and  denser  and  it  is  found  to  satisfy 
a  Rayleigh-Taylor  instability  criterion  at  -  40  nsec.  This  is  roughly  the 
time  when  the  rearward  accelerated  debris  has  swept-up  an  equivalent  target 
mass.  The  overall  results  further  suggest  that  the  radiating  line  cores 
should  be  optically  thick. 

In  the  next  seobion  we  discuss  the  code  in  detail  by  enumerating  the 
physical  processes  included  in  it.  In  section  III  we  describe  a  standard  run 
and  make  comparisons  with  other  runs  which  employ  ad  hoc  assumptions  for 
chemistry  and  radiative  line  transfer.  In  the  final  section  we  summarize  our 


main  conclusions  and  ma<e  ccmoarisons  with  the  exoeriment.  The  limitations  of 


he  present  code  are  also  presented  here. 
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II.  DISCUSSION  OF  THE  SIMULATION 


A.  Hydrodynamics. 

The  basic  hydrodynamic  equations  to  be  solved  are  the  continuity 
equation, 
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the  total  momentum  equation, 
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the  ion  internal  energy  equation, 
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and  the  electron  internal  energy  equation. 


3  t  ^ 
3t(Ve} 


7*VpeEe  +  peV#v 


V.q£ 


-Qie  + 


A  ♦  H 


(«) 


Here  p  is  the  total  mass  density,  pg  is  the  electron  mass  density,  v  is  the 

mass  velocity,  p^  and  pg  are  the  ion  and  electron  pressure,  qi  and  q^.  are  the 

ion  and  electron  thermal  conduction,  z.  and  e  are  the  ion  anc  electron 

i  e 

specific  internal  energies,  A  is  the  electron  energy  change  due  to  radiation 
and  ionization,  and  H.  is  the  electron  heating  due  to  absorption  of  the  laser 
energy  to  be  described  below  (eqn.[30]).  Qie  is  the  thermal  equilibration 
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term  from  Braginskii 
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where  m.  is  the  mean  ion  mass  and  t  is  the  electron  collision  time.  Q  ,  _ 
i  e  vis 

is  a  tensor  artificial  viscosity  added  to  smooth  out  shocks  and  it  has  the 
form  suggested  by  Mch^  including  an  artificial  heat  flux  near  shocks.  A 
tensor  product  involving  the  artificial  viscosity  is  denoted  by  a  colon.  We 
have  also  assumed  that  there  are  no  currents  so  that  the  ion  and  electron 
velocities  are  equal. 

We  will  limit  our  analysis  so  that  the  variables  depend  on  time  and  only 
one  spatial  coordinate  in  a  general  orthogonal  coordinate  system.  Let  this 
coordinate  be  51  and  the  scale  factors  be  h, ,  h,,  and  h^.  Then  with  v  -  ve1 , 
q  -  qe. ,  and  Q  -  Qe^e  one  finds  for  the  continuity,  momentum,  and  internal 
energy  equations. 
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Let  us  define  the  area  as 
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where  the  integration  limits  extend  over  the  allowed  range  in  and  Mow 
note  that  for  an  arbitrary  function  Xf^), 
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where  the  range  of  integration  over  £  is  a  to  b  and  D/Dt  is  the  time 
derivative  following  the  fluid.  Integrating  each  equation  ( 6 )— ( 1 0 )  over  a 
volume  element  and  using  the  above  relation  gives 
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Squaticns  (’l)-(iti)  are  differenced  on  a  Lagrar.gian  grid  where  p,  pg, 

e.t  s  ,  o.,  and  p  are  carried  at  the  ceil  centers  and  v  i3  carried  at  the 
i  e  •  i  re 

ceil  interfaces.  There  are  as  many  of  the  above  equations  as  there  are  cells 
in  the  whole  grid  and  the  limits  of  integration  for  the  continuity  and 
internal  energy  equations  are  the  cell  interfaces,  while  the  limits  are  the 
ceil  centers  for  the  momentum  equation.  The  differencing  is  second  order  in 
space  and  a  two  step  predictor-corrector  scheme  is  used  for  the  temporal 
advance . 


Identifying  the  Kinetic  energy  of 
as 


ceil  J+1/2  with  interfaces  at  5.  and 
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one  can  show  through  a  lengthy  algebraic  computation  that  the  differencing 
scheme  conserves  total  energy  even  though  the  internal  energy  is  transported. 
This  is  valid  as  long  as  the  velocities  at  the  boundary  interfaces  are  zero. 

Due  to  the  large  change  in  volume  during  the  experiment  from  the  target 
dimensions  (a  few  microns)  to  the  blast  wave  formation  (-1  cm),  it  is 
necessary  to  rezone  the  grid.  This  is  done  each  time  step  by  first 
determining  the  ceils  which  provide  the  mcst  stringent  Courant  condition.  An 
isolated  subset  of  these  cells  i3  chosen  and  the  interface  between  each 
selected  cell  and  its  neighbor  is  moved  in  suer,  a  way  as  to  relax  the  lour  ant 
condition.  The  amount  by  which  the  interface  is  moved  is  determined  by 
conserving  total  mass,  linear  momentum  and  energy. 

Finally,  we  mention  that  the  different  terms  in  the  above  equations  are 
handled  in  a  time  split  fashion.  The  hydrodynami os ,  including  the  artificial 
viscosity  are  cone  first,  then  the  thermal  conduction  and  equilibration  terms 
are  solved  implicitly  each  time  step.  With  these  terms  completed  the 
chemical  ’"“actions  are  next  computed  using  an  improved  AHI3"  modui.e,  and 
"adiatio.n  is  transported . 


finally  tne  r 


B.  Geometry. 


The  initial  configuration  of  the  experiment  is  a  thin,  planar  disc  target 
which  is  hit  on  one  side  with  a  laser  beam.  The  blowoff  from  the  ablation 
then  undergoes  a  divergent  expansion  into  the  background  gas.  Even  though 

the  code  is  one  dimensional,  i.e.,  the  variables  depend  on  only  one 
coordinate,  one  would  like  to  properly  treat  both  the  initial  configuration 

and  the  divergent  blowoff.  Clearly  a  Cartesian  geometry  can  handle  the 

planar  target  but  not  the  divergent  blowoff.  A  spherical  geometry,  as  usee 
7 

by  Stall ingwerf  in  the  code  MACH1  to  model  the  laser  target  interaction, 
does  account  for  the  divergent  expansion,  but  it  has  other  problems.  Given 
that  the  target  material  expands  in  a  cone  of  half-angle  e  ,  it  is  impossible 
for  the  section  of  the  sphere  representing  the  target  to  simultaneously  have 
the  same  area,  total  mass,  density  and  thickness  of  the  actual  target. 
Moreover,  the  most  limiting  factor  is  the  inability  of  a  spherical  geometry 
to  account  for  the  acceleration  and  divergence  of  the  rear  side  of  the  target 
after  tne  laser  radiation  has  been  absorbed.  In  3uch  a  geometry  the  rear 
side  implodes  toward  the  center  of  the  coordinate  system.  Shadowgraphs 
clearly  show  that  the  rear  side  blowoff  is  quite  different  from  that  of  the 
front  side  and  it  is  therefore  of  interest  to  model. 

A  spheroidal  oblate  coordinate  system  offers  the  advantage  of  both  a 
planar  initial  target  and  a  spherical,  divergent  expansion.  Such  a 
coordinate  system  can  be  generated  by  drawing  ellipses  in  a  plane 
perpendicular  to  the  target  face.  These  ellipses  converge  to  a  bar  at  tne 
target  face  and  tend  toward  a  circle  far  from  the  target.  The  orthogonal  set 
of  coordinate  lines  have  an  hourglass  shape  near  the  target  and  become 
spherical  rays  at  large  distances.  These  generating  lines  are  then  rotated 
about  an  axis  through  the  center  of  the  disc.  The  axis  represents  the  path 
of  the  laser  beam.  One  of  these  hourglass- shaped  generating  lines  is  chosen 


as  the  lateral  boundary  of  the  expanding  material.  Figure  1  is  a  sketch  of 
the  geometry  and  initial  configuration.  The  coordinate  transformation  to 
cartesian  coordinates  is  given  by 
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where  0  £  s  <  ®,  -w/2  £  tp  £  ir/2,  and  0  £  p  £  2t.  In  terms  of  the  discussion 

Q 

of  oblate  spheroidal  coordinates  given  by  Arfken  ,  s  -  a  sinh(u).  The  flow 
variables  are  taken  to  depend  on  the  distance  s  and  a  fixed  pQ  determines  the 


hourglass-shaped  lateral  boundary. 


From  early  observations  of  the 


experiment,  the  major  portion  of  the  target  material  expands  into  a  cone  of 

Q  Q 

half-angle  40  .  Thus  -pQ  -  40  ,  and  the  radius  of  the  target  disc  (r  )  that  is 
irradiated  by  the  laser  is  related  to  the  parameter  a  by 
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The  corresponding  scale  factors  are 
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From  these  the  elemental  area  and  volume  are  given  by 
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C.  Thermal  Conduction. 


The  electron  thermal  conduction  has  long  been  known  to  present  a 

q 

problem  in  laser-target  interactions  .  The  difficulty  lies  in  the  fact  that 
the  temperature  gradient  length  scale  -  |dir.(T  )/dx j  is  on  the  order  of 

the  mean  free  path  X  near  the  critical  surface,  while  the  classical  formula 
for  the  thermal  conductivity  assumes  x/LT  <<  1.  The  usual  approach  in  fluid 

codes  to  rectify  this  problem  has  been  to  limit  the  heat  flux  q  such  that 
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is  the  free 


where  «.  is  the  classical  or 

Ci. 

streaming  limit  for  thermal 
typically  0.03  <  a  <  0.1. 

achieved  by  using  instead 


3pitzer-Harm  conductivity,  and  qj«s 

heat  flux,  qfs  -  -VbVWV"2  with 
A  smooth  variation  in  the  heat  flux  can  be 
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with  <  »  <  ,  . 
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It  has  also  been  shown  by  Manheimer11  that  the  heat  flux  can  be  reduced 
by  turbulence  through  ion  acoustic  waves.  For  "T  ^  <<  T  ,  these  waves  are 

driven  unstable  by  the  return  current  of  an  anisotropic  electron  distribution 

1  p 

function.  The  criterion  for  instability  can  be  written  as 


k  T  1/2r  m  1/2 

q  >  F  -  4n.k  T  (-SL2)  ( 1 *3T . /T  )  /2(-^) 

2  e  3  e  m  ie  m. 

e  t  i 


(^)1/"[(l+3Ti/Te)1/2  -  (T./T,)172]  exp(-  -7;^-  ) 


1,(22) 


When  this  condition  is  satisfied  the  heat  flux  can  be  written  as 


1 1 


‘cl 


-7 - 77  , 

i  *  <./<  e 


(23) 


cl  an 


where  <  is  the  anomalous  thermal  conductivity  due  to  the  turbulence  3nd  is 


given  by 


with  f(Z)  a  function  of  the  mean  charge  state  Z,  XQe  is  the  Debye  length,  and 
<p  the  potential  fluctuations  due  to  the  turbulence.  The  instability 


saturates  by  ion  trapping  with  the  maximum  fluctuation 


e<*max  1 ,  1  -  J_i.2 

kBTe  ’  1.25  '  *  Te; 


(25) 


The  procedure  we  use  to  determine  instability  is  more  efficient  than  that 

described  ir.  Ref.  11  or  Ref.  12.  From  eons. (22)  and  (23)  we  first  determine 

the  sign  of  a  «  q  /F  -  where  qQl  -  icc^7T0  30(3  *  3s  de^*ned  in  eqn.(22). 

If  a  _<  0,  then  the  waves  are  stable  and  k  -  k^.  On  the  other  hand,  if  o  is 

positive  then  the  instability  is  operative.  Next  the  magnitude  of  the 

fluctuations  is  found  by  setting  \|»  -  min  (q  .  -  1,  ^  ),  where  ^  is  given 

max  ma x 

by  eqn.tZM  using  eqn.(25).  This  procedure  accounts  for  the  saturation  in 
the  growth  of  the  unstable  waves  and  gives  the  marginal  stability  condition 
otherwise. 

The  effect  of  saturation  of  the  heat  flux  in  steep  temperature  gradients 

must  still  be  included,  so  once  having  found  the  value  of  ie  , /<  we  use  tc  - 

cl  an 

/(I  +  <_,/<,_)  in  eon. (21).  In  the  running  of  the  code  it  was  found  that 
vi  vX  an 

the  heat  flux  saturated  near  the  critical  surface  but  the  ion  acoustic 
instability  was  stable  there  due  to  the  fast  thermal  equilibration  times.  At 
and  in  front  of  the  blast  wave  where  the  thermal  conduction  leads  to  a 
leading  electron  temperature  foot,  the  electron  heat  flux  was  limited  by  ion 
acoustic  turbulence  as  well  as  being  saturated. 

A  third  effect  related  to  the  heat  flux  is  its  non-local  nature  in 
extreme  temperature  gradients.  This  aspect  has  beer,  discussed  by  Luciani,  et 
al J  ^ ,  bUt  was  not  included  in  the  present  code.  Further  work  on  this 
problem  in  the  presence  of  large  density  changes  is  needed  before  their 
expression  for  the  heat  flux  can  be  used  ir.  the  present  simulation. 
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We  finally  mention  that  ion  heat  flux  is  included  in  the  code  and  the 
effect  of  saturation  is  treated  in  a  similar  manner  to  that  of  the  electron 
heat  flux.  We  did  not  consider  any  instability  which  could  limit  the  ion 
heat  flux  by  turbulence. 

D.  Equation  Of  State. 

The  equation  of  state  is  different  from  that  of  an  ideal  gas  law  only  in 
the  target  material  while  it  remains  a  solid  prior  to  vaporization.  Although 
the  energy  required  to  completely  vaporize  the  target  foil  is  very  small 
compared  to  the  total  laser  energy  in  the  experiment,  neglect  of  a  realistic 
equation  of  state  for  the  solid  target  leads  to  absurd  compressions  in  the 
rear  side  of  the  target.  We  do  not  use  a  detailed  model  or  table  for  the 
solid  equation  of  state,  for  this  would  require  spending  an  excessive  amount 
of  computer  time  on  transient  phenomena.  Rather  we  employ  a  simple  analytic 
model  which  adequately  accounts  for  the  dominant  physical  effect,  namely  the 

rapid  rise  in  the  sound  speed  in  a  compressed  solid. 

.  1 U 

ihe  totax  pressure  can  be  written  as 


p  •  pc  *  pn  *  pe  •  (26) 

where  subscript  c  refers  to  the  zero  temperature  isotherm  or  ccld,  solid 
component,  subcript  n  to  the  nucleon  component  (neutrals  or  ions),  and 
subscript  e  to  the  electron  component.  Let  be  the  density  of  the  solid 

at  zero  temperature  and  n  »  p/pQ0>  then  the  internal  energy  and  pressure  are 
related  by 


P„n  2 


dn  . 


(27) 


Based  on  formulas  (3.4)  and  (3.20)  in  Ref. lit  we  employ  the  following 
expressions  for  the  cold  pressure: 


c32(n5/3-  1) 


,  n  >  1 


pc  “ 


C,n2/3C 


c_(i-n"1/3) 


-  e 


c60-n 


■1/3. 


]  ,  n  <  1 


(28) 


-1  1  3 

where  for  aluminum  e  -  ergs/cm,  C_  -  7.22,  C,  -  3.Uf  and  p 

32  00  00 

■5  , 

2.71  g/cm  .  The  value  of  was  taken  to  be  much  smaller  than  c.l*1C 

■3  *73 

ergs/cnr  of  Ref.  It,  namely,  C„  -  1*10  ergs/cm  .  As  can  be  seen  in  the  P- 

V-T  diagram  for  aluminum  in  Fig.  D-8  of  Ref.  14,  there  is  a  deep  but  narrow 

well  of  negative  pressure  for  r\  slightly  less  than  1.  This  feature 

represents  the  intramolecular  forces  which  bind  a  solid  together.  A  large 

negative  pressure  can  lead  to  numerical  instabilities  unless  exceedingly 

small  time  steps  are  used.  Since  the  phase  of  negative  pressures  is 

transitory  and  the  internal  energy  associated  with  the  well  is  small  we  can 

reduce  the  depth  of  the  pressure  well  without  seriously  altering  the  physics. 

Further,  as  the  target  rarefies  due  to  heating  the  thermal  pressure  terms  in 

eqn.(26)  readily  dominate.  From  eqns.(27)  and  (28)  and  the  coefficients  C^, 

Cj.,  Cg  from  Ref.  14,  one  finds  that  the  zero  temperature  sublimation  energy 

+4 

for  the  target  geometry  of  Fig.  1  is  3.3*10  ergs,  which  is  much  less  than 
the  23  Joules  deposited  by  the  laser. 

The  steep  rise  of  the  pressure  for  compressions  does  play  a  significant 
role  in  the  propagation  of  shocks  through  the  solid  target.  In  agreement 
with  the  results  of  Schmaiz  and  Meyer-ter-Vehn1 ^  who  used  a  detailed  equation 
of  state,  we  also  found  a  series  of  shocks  in  the  target  during  the  early 
stages  of  disassembly.  These  shocks  display  a  small  density  jump  but  large 
pressure  changes. 


r'/’-S 
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The  other  terms  in  eqn.(26)  are  treated  according  to  the  ideal  gas  law, 


except  initially  when  the  temperature  is  <  0.1  ev.  The  P-V-T  diagram 

indicates  a  constant,  low  pressure  for  low  temperatures  and  n  <  1.  In  these 
regions  the  pressure  in  the  solid  is  taken  to  be  that  of  the  initial  ambient 
gas.  This  is  reasonable  since  the  foil  and  background  gas  are  in  equilibrium 
before  the  laser  is  turned  on,  while  the  temperature  rapidly  increases 
af  ter  wards . 


S.  Radiation. 

We  begin  with  a  discussion  of  the  initial  radiation  in  the  problem,  that 
cf  the  incident  laser.  The  total  energy  in  the  laser  pulse  was  set  at  23 
Joules  and  a  Gaussian  shape  of  4  nsec  full-width  half-maximum  was  adopted. 
The  laser  wavelength  was  1.06  uo,  modeling  the  output  from  a  neodymium  glass 
laser.  The  peak  intensity  of  the  pulse  was  IkIO*1^  W/cm'  and  occured  5  nsec 
after  the  3tart  of  the  calculation.  The  pulse  was  cut  off  at  10  nsec.  The 
beam  energy  was  deposited  using  the  1-0  radiative  transfer  equation  with 
inverse  bremsstrahlung  absorption: 


< 


IB 


(■i)  / >4).  )  ' 

oe  L 


<.U)  '  -W*  ' 

pe  u 


.1/2 


(29' 


where  ui  is  the  electron  plasma  frequency.  is  the  electron- ion  collision 
time,  and  m  is  the  laser  frequency.  Since  the  radius  of  t re  beam  was  that 
of  the  target  di3c  (r^),  the  heating  rate  per  'unit  volume  for  the  electrons 


in  cell  J  is 


w. 


-  -  - 


d  *J+l 


(1  -  e  ) 


Vol ( Jj 


where  I .  is  the  mean  intensity  entering  cell  J  of  thickness  A . .  The 
J  '  J 

fractional  ionization  at  the  front  edge  of  the  target  is  artificially 
increased  at  the  initial  instant  of  the  simulation  to  provide  free  electrons. 
This  ensures  the  absorption  of  the  laser  beam  since  inverse  bremsstrahlung 
depends  on  the  presence  free  electrons.  For  simplicity,  we  neglect 
reflection  of  the  laser  radiation  and  deposit  any  remaining  incident 
radiation  not  absorbed  before  the  critical  surface  into  that  cell  containing 
the  critical  surface.  The  critical  surface  is  defined  as  the  position  where 
the  denominator  of  eqn.(29)  vanishes. 

The  radiation  from  the  heated  gas  is  treated  in  several  ways.  First, 

bremsstrahlung  and  recombination  radiation  are  calculated  for  both  the  target 

aluminum  material  and  the  ambient  nitrogen  gas.  The  recombination  radiation 

is  assumed  to  have  a  Planck-like  spectrum  with  an  effective  temperature  T„ad. 

For  cell  J  this  quantity  is  determined  by  equating  the  energy  density  of  a 

blackbody  at  T  ,  with  the  energy  emitted  by  recombination  radiation  in  cell 
raa 

J.  Radiation  emitted  by  the  nitrogen  gas  is  treated  as  optically  thin,  i.e., 
an  energy  less.  However,  radiation  from  the  aluminum  is  transported  through 
the  entire  grid.  The  radiation  transport  is  calculated  by  the  standard  ray¬ 
tracing  technique  in  the  forward  and  rearward  direction.  As  shown  in  the 
detailed  paper  by  Duston  et  al.1^,  radiation  is  a  major  source  of  heating  for 
the  rear  side  of  the  target  foil.  Hence,  since  one  of  our  objectives  is  to 
simulate  the  rearward  blowoff,  it  is  imperative  to  correctly  calculate  the 
basic  energy  transport  through  the  entire  target.  The  radiated  energy  is 
divided  into  twelve  frequency  bins  which  are  chosen  to  resolve  the  inner- 
shell  absorption  edges  of  aluminum.  The  absorption  cross  sections  are  from 
Fig.  2  of  Ref.  16.  Any  radiation  that  passes  through  the  aluminum  gas  is 


tv 
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also  transported  through  the  ambient  nitrogen  gas.  This  results  in  the 

dissociation  of  background  M^.  The  absorption  cross  sections  for  this 
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process  are  based  on  those  used  in  Hyman  et  al. 

Second,  excitation  line  radiation  is  modelled  in  two  extreme  conditions: 
completely  optically  thick  and  optically  thin.  The  optically  thick  case 
simply  means  that  any  emitted  line  from  a  cell  J  is  absorbed  in  that  same 
cell.  No  calculation  is  needed  here.  For  the  optically  thin  case  we 
calculate  the  emission  from  13  resonance  lines  of  aluminum  and  14  resonance 
lines  of  nitrogen.  All  ionization  stages  of  both  species  are  covered.  Cur 
interest  in  the  line  radiation  stems  from  the  possibility  of  cooling  in  the 
cavity  behind  the  blast  wave.  We  did  not  consider  the  X-lines  which 
originate  during  the  early  phase  of  target  vaporization.  In  calculating  the 
emission  from  cur  chosen  set  of  lines,  the  excited  state  populations  are 
determined  by  balancing  the  collisional  excitation  rate  against  the  de¬ 
excitation  rate  from  the  excited  state,  the  ionization  rate  to  the  next 

higher  ionization  stage,  and  the  radiative  decay  rate.  Optical  pumping  as 
well  as  photoionization  of  the  resonant  state  due  to  recombination  radiation 
are  not  evaluated.  These  features  ensure  that  the  line  radiation  is  not 
overestimated  in  the  absence  of  reabsorption. 


7.  Chemistry. 

Much  effort  was  spent  in  developing  a  chemistry  module  which  can  handle 
the  extreme  range  of  chemical  conditions  in  the  experiment.  For  simplicity, 
we  do  not  follow  any  excited  state  populations  but  assume  all  ions  {ana 
neutrals)  are  in  the  ground  state.  Essentially  we  divided  the  chemical 
evolution  of  tr.e  aiummun  material  into  three  schemes.  First,  during  a  brief 


initial  stage  the  aluminum  undergoes  a  phase  transition 


from  a  solid  to  a 


high  density  gas  at  <  1  ev.  Here  we  use  an  analytic  procedure  to  determine 
the  temperature  and  fractional  ionization  from  the  internal  energy.  The 
analytic  solution  is  based  on  fits  to  numerical  calculations  of  an 
equilibrium  aluminum  gas  with  only  neutral  and  singly  ionized  species.  As 
the  target  continues  to  absorb  the  laser  energy,  the  fractional  ionization 
and  temperature  drastically  increase.  Once  the  fractional  ionization 
increases  above  0.8,  an  iteration  on  the  temperature  and  ionization  level  for 
a  given  internal  energy  is  employed,  again  for  equilibrium  conditions.  In 
this  second  scheme  for  the  chemical  evolution  all  the  ground  state  ionization 
levels  are  included.  Finally,  when  the  electron  density  in  any  cell  falls 
below  5. *10  cm  J,  the  chemistry  in  that  ceil  is  integrated  in  time  using 
reaction  rates,  i.e.,  we  allow  the  population  of  the  ground  state  levels  to 
fall  out  of  equilibrium.  For  this  third  scheme,  the  ARIS^  module  is  used  to 
solve  the  corresponding  set  of  stiff  differential  equations.  For  the 
nitrogen  gas,  the  chemical  evolution  is  always  followed  using  reaction  rates. 

The  numerical  details  of  this  evolution  from  low  temperature,  dense 
equilibrium  chemistry  to  high  temperature,  rarefied  reaction  rate  chemistry 
is  too  complicated  to  be  reported  here.  We  do  mention  that  the  transition 
from  one  chemistry  scheme  to  the  next  must  be  done  smoothly.  This  requires 
the  inclusion  of  many  reaction  processes  to  ensure  that  the  limiting 
conditions  lead  to  equilibrium  results.  For  instance,  at  high  density  and 
lew  temperature  we  reproduce  the  simple  Saha  equilibrium  by  using  detailed 
balance  between  collisional  ionization  and  3~body  recombination.  Such 
equilibrium  is  valid  only  for  a  high  density  gas  where  the  dominant 
recombination  process  is  3-body.  In  addition  to  3~t>ody  recombinations  and 
electron  collisional  ionization,  we  also  include  radiative  recombination  and 
dielectrcnic  recombination  processes  so  that  the  lower  density,  higher 
temperature  regimes  are  modeled  correctly.  A  high  density  correction  for 
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dielectrcnic  recombination  rates  similar  to  that  used  by  Jordan  is 
employed. 

As  an  aside  on  computational  physics,  the  use  of  a  sequence  of  numerical 
schemes  for  the  chemistry  is  most  efficiently  and  correctly  employed  in  a 
Lagrangian  hydrodynamic  code.  In  this  approach  the  decision  to  switch  from 
one  scheme  to  the  next,  vis.,  analytic  ->  equilibrium  ->  time  dependent,  can 
be  made  in  each  individual  cell  independent  of  all  the  others  since  each  cell 
carries  the  same  material  with  it.  Unfortunately,  in  a  purely  Lagrangian 
code  the  effects  of  charge  exchange  cannot  be  investigated,  for  the  aluminum 
debris  cannot  overrun  the  nitrogen  background.  We  must  thereby  limit  the 
simulation  results  to  the  high  pressure  regime  where  strong  coupling  between 
the  expanding  debri3  and  background  occurs. 


SAMPLE  HUNS 


In  this  section  we  will  present  detailed  graphical  results  from  the 
numerical  simulations.  Four  distinct  models  were  studied.  The  differences  in 
the  models  lies  in  the  treatment  of  the  chemistry  and  line  radiation.  In 
model  A,  which  we  refer  to  as  the  standard  model,  time  dependent  chemistry  is 
employed  and  all  line  radiation  is  assumed  to  be  optically  thick,  i.e.,  no 
transport  of  line  radiation  is  considered.  In  model  5,  equilibrium  chemistry 
is  used  and  again  the  line  radiation  is  optically  thick.  For  models  C  and  D 
the  line  radiation  is  taken  to  be  optically  thin,  but  model  C  has  time 
dependent  chemistry  and  model  D  has  equilibrium  chemistry.  Recall  from 
section  II. F.  that  equilibrium  chemistry  is  net  necessarily  LTE  but  includes 
the  effects  of  radiative  and  dielectrcnic  recombination.  The  objective  of 
studying  these  four  models  is:  (i),  to  demonstrate  the  importance  of  non¬ 
equilibrium  chemistry,  and  (ii),  to  investigate  whether  line  radiation  can 
realistically  cool  the  cavity  inside  the  blast  wave.  Since  we  do  not  perform 
a  detailed  calculation  for  the  transfer  of  line  radiation,  the  two  extremes 
of  optically  thick  and  thin  enable  one  to  obtain  only  qualitative 
conclusions . 


A.  Standard  Model. 


Figure  2  3hows  the  position  of  the  blast  wave  for  both  the  front  side 
and  the  back  side  as  a  function  of  time  for  model  A,  the  standard  model.  The 
time  in  this  figure  is  measured  from  the  instant  of  peak  intensity  of  the 
laser  pulse,  which  is  what  the  experimentalists  use  as  their  zero  time.  In 
the  code  the  laser  pulse  is  initiated  5  nsec  before  the  peak.  The  solid  line 
in  Fig.  2  is  given  by 


[t(nsec)] 


R(om) 


0.123  [ 


£(  Joules )  -|0.2 
p(Torr) 


0.4 


(31  ) 


for  E-23  Jcules  and  p  -  2.5  Terr.  This  equation  i3  from  Ref.  3  anci 
describes  the  observed  temporal  evolution  of  the  front  side  of  the  blast 
wave.  As  noted  in  that  reference,  the  equation  for  R(t)  has  the  sane 

temporal  and  parameter  dependence  as  that  for  a  Taylor-Sedov  blast  wave.  We 
note  the  excellent  agreement  in  Fig.  2  between  the  numerical  simulation  and 
the  experimental  data  for  the  position  of  the  front  side  blast  wave.  The 

blast  wave  on  the  back  side  displays  a  turbulent  structure  in  the 

shadowgraphs  and  no  general  formula  for  it3  position  in  time  has  been 

developed.  This  turbulent  structure  will  be  discussed  below  in  light  of  the 
results  of  the  simulation.  The  equal  mass  radii  (R_M)  for  both  the  front  and 
back  side  are  also  noted  in  Fig.  2.  Specifically,  0.15*10~6  g  is  ablated 
toward  the  front  side  while  the  remaining  0.55*1 g  is  accelerated 

rearward.  As  demonstrated  by  the  simulation,  the  Taylor-Sedov  blast  wave 
solution  becomes  applicable  for  the  front  side  shock  location  only  after 

several  equal  mass  radii  rave  been  encountered. 

In  Fig.  3  the  energy  budget  for  model  A  is  presented.  ,'E.  dt  is  the 

U 

accumulated  energy  deposited  by  the  laser  over  time.  It  attains  23  Joules  at 
*0  nsec  after  the  peak  of  the  laser  pulse.  £..  is  the  instantaneous  kinetic 
energy  of  directed  motion  while  is  the  instantaneous  thermal  (random) 

energy  of  all  particles.  The  chemical  energy,  i.e.,  tr.e  energy  looxed  up  in 
the  ionization  of  aluminum  and  nitrogen  at  a  given  time  is  denoted  by 
Finally,  /E^dt  is  the  accumulated  energy  lost  from  the  system  during  the  run 
due  to  recombination  and  bremsstranlung  radiation.  The  values  in  Fig,  3 
include  both  sides  of  t.ne  expanding  target.  The  relation 


s.,(c)  *  E„(t;  -  £„'t ;  - 


£„dt 


expressing  conservation  of  total  energy,  is  satified  to  within  3$  tnroughout 


the  run.  This  resuit  is  also  typical  for  the  other  runs  dicussec  below. 
Nearly  all  of  the  energy  that  is  missing  during  a  simulation  run  occured  in 
the  time  dependent  chemistry  routine. 

Note  in  Fig.  3  that  initially  the  kinetic  energy  is  dominant,  but  after 
-  +20  nsec  the  fraction  in  thermal  energy  is  largest  and  stays  at  roughly 
35$.  For  an  adiabatic  Taylor-Sedov  blast  wave  with  Y  -  1.1  (see  ref.  3)  90$ 
of  the  total  energy  is  thermal  with  the  rest  being  kinetic.  This  is  clearly 
not  the  case  here  and  points  to  the  first  major  difference  between  the 
simulation  of  the  experiment  and  the  ideal  blast  wave  theory.  This 
difference  is  due  to  the  significant  amount  of  energy  locked  up  in  ionization 
potential  and  the  smaller  radiated  energy.  These  two  energy  sinks,  which  are 
not  included  in  the  ideal  Taylor-Sedov  model,  account  for  roughly  45$  of  the 

total  energy  in  the  system.  The  high  fraction  of  kinetic  energy  at  the  end 

of  the  run,  20$  instead  of  10$,  arises  from  the  fact  that  the  rear  side 
blowoff  is  dense  and  relatively  cool. 

Figures  4a  through  4g  show  snapshot  profiles  of  different  quantities 
during  the  standard  model  simulation.  The  figure  captions  fully  describe  the 
nomenclature  and  symbols  used  to  present  the  numerical  data.  Here  we  will 
point  out  a  few  salient  features. 

Note  first  of  all  that  the  time  shown  on  the  graphs  refers  to  the 

simulation  time  which  is  zero  when  the  laser  turns  on.  To  transform  to  the 

time  from  the  peak  intensity  of  the  laser  pulse  subtract  5  nsec.  Figure  4a 

displays  the  situation  at  the  initial  instant  of  the  simulation;  the  aluminum 

is  at  solid  density  of  6*10  atoms/cnr  and  the  ambient  nitrogen  is  at 
+  17  3 

3.5*10  molecules/ cm'.  3oth  materials  are  predominantly  neutral,  but  a 
thin  partially  ionized  region  at  the  front  of  the  target  is  imposed  to  ensure 
that  the  laser  beam  is  absorbed  by  inverse  bremsstrahlung.  Note  that  the 
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abscissa  is  in  Lagrangian  cell  units  in  order  to  expand  the  target  region  j.  or 
clarity. 

Figures  lb,  4c ,  and  Id  display  the  situation  shortly  after  the  peak  of 
the  laser  pulse  again  in  Lagrangian  cell  units.  The  critical  surface  occurs 
near  the  peak  of  the  electron  density  and  temperature  (see  Fig.  lb).  One  can 

dearly  see,  from  the  plot  of  the  artificial  viscous  pressure  Q,  a  series  of 
shocks  traversing  the  aluminum  target.  This  feature  agrees  with  the  detailed 
equation  of  state  study  of  Ref.  15.  The  heating  cf  the  back  side  of  the 
target  is  due  to  these  shocks  as  well  as  to  the  absorption  of  radiation 
emitted  from  a  region  near  the  critical  surface.  The  graph  of  the  mean 
charge  state  in  Fig.  ic,  n  /(n.  ♦  n  ),  shows  that  highly  ionized  aluminum, 
namely  AZ+10,and  Al+11,  is  created  early  on  in  the  laser  target 

interaction.  We  find  that  during  the  subsequent  stages  cf  forward  blowoff 
ana  cavity  formation  the  aluminum  does  not  have  enough  time  to  recombine  and 
remains  in  these  highly  ionized  states.  Since  the  cavity  temperature  drops 
below  TOO  ev,  this  non-equilibrium  ionization  structure  points  cut  the 
importance  of  using  a  non-LTE,  time  dependent  chemistry  in  the  simulation. 
The  spectrum  in  Fig.  id  only  includes  the  radiation  from  recombination  and 
bremsstra.niung;  it  displays  some  of  the  qualitative  features  of  Figure  11  in 
Ref. 16  but  clearly  the  lack  of  a  detailed  treatment  of  the  recombination 
spectrum  and  lines  represents  a  limitation  of  our  present  model. 

Figure  ue  is  a  snapshot  cf  the  flow  variables  at  the  time  when  the 
forward  moving  ol owcr'f  has  swept  up  one  equal  mass  of  ambient  material.  The 
abscissa  is  r.ow  given  in  physical  space  since  the  target  has  expanded 
sufficiently  to  ce  resolved,  ""he  formation  cf  a  reverse  shook  at  this  time 
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Figures  4f  and  4g  show  the  long  time  development  cf  the  system  out  till 

the  forward  moving  blast  wave  reaches  the  boundary  of  the  calculations!  grid. 

From  the  density  plot  in  Fig.  4f  one  finds  that  the  blast  wave  thickness,  as 

measured  by  the  width  of  the  peak  in  the  electron  density  profile,  is  roughly 

twice  as  thick  as  the  experimentally  observed  value  of  3$  of  the  blast  wave 

radius^.  However,  the  peak  electron  density  in  the  shock  of  4*10+1^  cm-^  is 

1 3 

typical  of  the  reported  measurements  from  several  interf erograms  . 
Moreover,  the  mean  electron  density  in  the  cavity  of  a  few  times  10  cm  J 
is  in  good  agreement  with  the  recent,  though  tentative,  measured  value  of 

+  t  7  -■a  iq 

2*10  cm  J  at  60  nsec  after  the  peak  of  the  laser  pulse  .  The 

surprisingly  high  value  for  the  electron  density  reflects  the  high  degree  of 

ionization  of  the  aluminum  debris  in  the  cavity;  note  that  the  ion  density  in 

the  cavity  is  much  smaller.  From  the  viscous  pressure  plot  in  Fig.  4f  the 

reverse  shock  on  the  front  side  has  moved  back  toward  the  original  target 

position  and  the  formation  of  a  second  reverse  shock  in  the  backward  moving 

aluminum  is  also  evident.  3y  the  time  cf  Fig.  4g  the  cavity  between  the 

forward  and  backward  moving  blast  waves  is  nearly  isothermal  with  Tg  -  60  ev. 

This  is  the  second  major  difference  of  the  simulation  results  from  the  ideal, 

adiabatic  Taylor-3edov  blast  wave  wherein  the  temperature  rises  as  one  moves 

inward  from  the  blast  wave.  Actually  the  run  of  the  velocity  and  ion  density 

on  the  front  side  looks  quite  like  the  solution  for  an  isothermal  Tavlor- 

Sedov  blast  wave  as  discussed  in  Ref.  20.  The  isothermality  is  mainly  due  to 

electron  thermal  conduction.  We  note  for  reference  that  an  isothermal  blast 

0  4 

wave  also  evolves  according  to  a  t  ’  law. 

3.  Comparisons  With  Other  Models. 

Figures  5a,  5b,  and  5c  show  snapshots  of  the  flow  variables  for  models 


3,  C,  and  D,  respectively.  Each  of  these  graphs  is  at  a  time  close  to  50 


nsec  in  simulation  time  and  can  be  compared  with  rig.  4f.  The  distinction 
between  these  models  and  the  standard  one  (model  A)  is  noted  in  Table  I  which 
contains  some  of  the  gross  characteristics  for  each  model.  The  figures  and 
table  will  be  discussed  together.  The  main  differences  between  the  time- 
dependent  (model  A)  and  equilibrium  chemistry  (model  3)  are  the  lower 
electron  temperature  in  the  cavity  (65  ev  for  model  A  vs.  90  ev  for  model  3) 
and  the  higher  mean  ionization  level  in  the  cavity  (9.7  for  model  A  vs.  7.5 
for  model  3).  These  differences  arise  because  the  equilibrium  model  allows 


sharing  between  ionization  and  thermal  energy  with  no  regard  for  time  scales. 
When  we  assume  that  the  radiation  from  the  resonance  lines  is  optically  thin, 
as  in  models  C  and  D,  we  find  a  significant  increase  in  the  percentage  of 
radiation  lost  from  the  system.  The  values  of  61  %  and  38$  are  much  higher 
than  the  experimentally  estimated  25$  radiation  energy  loss  when  the  MRC  gas 
mixture  was  used  for  the  ambient  material.  Since  the  MRC  mixture  is  a  much 
more  efficient  radiator  than  a  pure  nitrogen  gas  we  surmise  that  the 
resonance  lines  cannot  be  optically  thin.  Simple  estimates  for  the  optical 
depth  of  a  Doppler  broadened  line  also  show  that  the  optical  depth  is  greater 
than  one  for  the  cavity  conditions.  Furthermore,  Figs.  5b  and  5c  indicate 


too  thin  a  blast  wave  and  too  high  a  peak  electron  density,  when  compared 


against  the  observations.  This  is  especially  true  for  model  Z. 
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TABLE  I 


COMPARISON  OF  DIFFERENT  MODELS 
■  23  Joules,  pQ  *  2.5  Torr ,  t  -  50  nsec) 

A  3  C  D 

Time  dep.  Equil.  Time  dep.  Equil . 

chem.  chem.  chem.  chem 

No  line  No  line  Line  rad.  Line  rad 


.  Further  Results. 


The  experimentally  observed  floculli  on  the  backside  blowoff  is  very 
distinct  from  the  smooth  front  side.  Within  the  confines  of  a  one 

dimensional  calculation  it  is  impossible  to  model  such  turbulence;  however, 
we  find  that  the  conditions  for  a  Rayleigh-Taylor  instability  are  satisfied 
in  the  back  side  blast  wave.  It  is  possible  that  such  an  instability  leads 
to  the  observed  structure  and  we  present  the  numerical  results  of  a  linear 
analysis  below. 

Specifically,  the  Rayleigh-Taylor  instability  arises  from  the 

deceleration  of  the  aluminum  debris  by  the  less  dense  ambient  nitrogen.  This 
deceleration  occurs  when  the  rearward  moving  debris  has  swept  up  about  an 
equal  mass  of  the  ambient  material.  The  growth  rate  fcr  the  decelerating 
Rayleigh-Taylor  instability  is  given  by 

»  -  C  Vin(p)  ]1/2  ,  (32) 

where  u  is  the  growth  rate.  This  expression  i3  calculated  at  each  interface 
between  cells  and  the  expression  for  the  temporal  velocity  derivative 
includes  the  pressure  gradients  and  the  artificial  viscosity  terms  in 
eqn.(12).  Whenever  <u  is  real  the  interface  i3  unstable.  The  resultant 
growth  rates  (to  real)  for  four  different  times  in  the  standard  model  are 

presented  in  Figure  6.  At  an  early  phase  in  the  simulation  the  front  side 

contact  discontinuity,  located  at  the  right  vertical  line  near  cell  ,  is 
seen  to  be  unstable.  However  the  duration  of  the  conditions  for  a  positive 
growth  rate  last  only  a  couple  of  e-folding  times,  and  by  55  nsec  the  front 
contact  discontinuity  is  essentially  stable.  (Isolated  crosses  are  not 
considered  to  be  indicative  of  consistent  conditions  for  Rayleigh-Taylor 
growth  but  rather  transient  phenomena.)  On  the  other  hand,  from  roughly  20 
nsec  to  at  least  55  .nsec  the  nearside  contact  discontinuity,  at  the  left 


vertical  line  near  cell  #20,  is  Ray! ei gh-T ayl or  unstable.  This  temporal 
duration  is  sufficient  for  several  e-folding  times.  As  can  be  seen  from  Fig. 
2,  this  time  span  centers  around  the  time  of  one  equal  mass  radius  for  the 
back  side  blowoff.  We  suggest  that  the  flocculi  observed  in  the  shadowgraphs 
at  -  50  nsec  and  beyond  are  due  to  this  Rayleigh-Taylor  deceleration 

instability.  The  continuation  of  the  turbulence  into  late  times  emerges  from 
the  sufficient  number  of  e-folding  times  for  this  instablility  on  the  rear 
side.  At  very  early  times,  say  10  to  20  nsec,  the  front  side  may  also 
exhibit  some  structure,  but  the  limited  number  of  e-folding  times  for  this 
region  would  not  lead  to  non-linear  effects. 

For  completeness,  we  display  the  local  Mach  number  at  four  different 

times  for  the  standard  model  in  Figure  7.  The  results  are  related  to  the 

21 

proposed  use  of  small  obstacles  as  a  probe  of  the  cavity  dynamics  .  Beyond 
about  20  nsec  the  major  portion  of  the  cavity  between  the  forward  and 
backward  moving  blast  waves  is  subsonic.  This  result  predicts  that  a  bow 
shock  formed  about  a  small  obstacle  at  the  time  it  is  overrun  by  the  forward 
moving  blast  wave  will  dissipate  into  a  sound  wave  as  the  expansion 
continues. 


IV.  SUMMARY 

We  have  presented  the  details  of  a  one  dimensional  numerical  simulation 
code  intended  to  model  the  NHL  laser/HANE  experiment  in  the  high  pressure 
regime  (p  >  1  Torr).  Although  the  code  is  one  dimensional  and  cannot  model 
the  development  of  instabilities  and  the  self-generated  magnetic  fields,  a 
large  number  of  physical  processes  have  been  modeled  in  order  to  study  the 
gross  physical  characteristics  of  the  experiment.  These  features  include  an 
oblate  spheroidal  geometry  to  account  for  the  rearward  expansion,  electron 
and  ion  thermal  conduction  with  saturation  including  the  limiting  effects  of 
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the  ion-acoustic  instability  on  electron  thermal  conduction,  an  equation  of 
state  representative  of  solid  aluminum  during  the  target  dissaseabiy , 
radiation  transfer  for  the  continuum  and  bremsstrahiung  radiation  originating 
in  the  aluminum  material,  and  an  intricate  chemistry  which  spans  the  range 
from  LTZ  tc  highly  non-equilibrium  conditions. 

Other  numerical  codes  have  studied  the  laser-target  interaction  problem 

in  2-0  with  accurate  radiation  transport  and/or  seif- generated  magnetic 
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fields,  such  as  that  of  Colombant,  et  al.  ,  Brackbiii  and  Goldman  ,  and 
Per t  .  These  cedes  are  constructed  to  simulate  the  first  few  hundred 
picoseconds  of  the  interaction,  before  any  hydrodynamic  reaction  of  the 
background  plasma.  The  emphasis  of  the  present  laser/HANE  experiment  i3 
primarily  on  the  coupling  of  the  debris  to  the  background  gas  which  takes 
place  over  a  timescale  of  tens  of  nanoseconds.  The  code  we  have  described 
concentrates  on  the  latter,  longer  timescale  phenomena,  and  hence  is 
significantly  different  from  the  previous  codes.  For  instance,  the  use  of  a 
Lagrangian  grid  or  adaptive  rezoning  i3  necessary  to  follow  the  great  change 
in  length  scales  from  the  target  size  to  one  centimeter:  the  above  referenced 
codes  employed  a  fixed  Sulerian  grid  for  the  hydrodynami os .  The  extension  of 
the  present  numerical  techniques  to  a  2-0  code  is  clearly  necessary  to 
improve  the  radiation  transport,  to  study  instabiiit.es,  and  to  include  the 
effeot3  cf  the  seif-generated  magnetic  field  on  thermal  conduction.  The 
present  mode!  .s  an  initial  3tep  toward  a  mere  complex  2~0  one  and  the 
results,  listed  below,  point  out  specific  features  which  a  more  advenced  code 


should  further  study. 

The  primary  results  f-ora  the  present  numerical  simulation  are  as 
follows . 


fi;  For  the  star.card  mode.  [A)  the  cavity  behind 
blast  wave  13  hot  'T  -  60  ev)  and  dense  (r.  -  2*10 

A  0 

temperature  in  the  cavity  is  accut  the  same  as  Ta  out  the 
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aluminum  part  of  the  cavity  is  a  factor  cf  10  smaller  (see  Fig.  *Jf). 

Overall,  the  dynamic  structure  of  the  forward  blowoff  region  resembles  an 

isothermal  Taylor-Sedov  blast  wave,  albeit  with  a  significant  amount  of 

energy  tied  up  in  chemical  potential  and  radiation  losses  (see  Fig.  3).  This 

<*»  ^ 

result  is  quite  different  from  the  LIE  model  presented  by  R.  Stellingwerf £ 
wherein  the  cavity  is  much  colder  (T  -  few  ev)  and  rarefied  (n.  -  few  *  101® 

-■3 

cm  J).  Our  predicted  electron  density  for  the  cavity  agrees  with  a  tentative 

*  Q 

measurement  by  J.  Stamper  .  The  electron  densities  in  the  blast  wave  also 
agree  with  the  measurements1®. 

(ii)  Due  to  the  initial  rapid  heating  and  the  slow  recombination  rates, 

the  material  in  the  cavity  is  highly  overionized  for  the  local  electron 

temperature.  In  particular,  by  50  nsec  the  aluminum  is  mostly  ,  A£,+1^, 

and  Al  ,  and  the  nitrogen  is  mostly  N  (  see  Table  I).  The  latter  has 
been  observed  in  an  integrated  spectrum.  and  it  should  be  possible  to  search 
for  the  presence  cf  AH.*10  using  the  same  technique.  Such  an  observation 
would  provide  a  direct  verification  of  our  use  of  a  time- dependent  chemistry. 
Due  to  the  limitation  cf  a  Lagrangian  code,  we  did  not  include  charge 
exchange  as  one  of  the  chemical  reactions.  If  the  blast  wave  becomes  leaky 
for  lower  ambient  pressures  (pQ  <  1  Torr),  seme  of  these  highly  ionized 

particles  could  penetrate  ahead  of  the  blast  wave.  To  model  charge  exchange 

and  leaky  blast  waves  a  two  ion  fluid  code  or  a  particle  code  would  be 

needed. 

(iii)  The  results  of  our  models  which  assume  optically  thin  lines 
(models  C  and  D)  indicate  an  unrealistically  large  amount  of  radiation  loss 
(see  Table  I).  This  suggests  that  the  resonance  lines  from  the  blast  wave 
and  cavity  are  not  optically  thin,  but  an  improved  radiation  model  with  line 
transport  is  required  to  determine  how  thick  the  lines  are. 

(iv)  The  deceleration  Rayleigh-Taylor  instability  is  found  to  be 
satisfied  at  the  rear  contact  discontinuity  at  the  time  when  an  equal  mass  is 


swept  up  by  the  rearward  moving  blast  wave  (see  Fig.  6).  We  suggest  that 
thi3  instability  may  be  responsible  for  the  flocculi  observed  on  the  rear 


side.  Again,  a  two  dimensional  code  i3  required  to  follow  the  non-linear 
development  of  the  unstable  regions.  On  the  front  side  the  instability  is 
also  active,  but  only  for  a  brief  time  near  7-20  nsec.  This  span  of  time 
allows  only  a  few  e-folding  times.  It  would  be  of  interest  to  take 
shadowgraph  pictures  at  these  early  times  to  discern  if  the  blast  wave  is 
initially  unstable  but  later  heals  itself  to  the  smooth  form  seen  in 
shadowgraphs  at  times  >  50  nsec.  It  should  be  noted  that  thi3  instability 
appears  to  have  little  to  do  with  the  formation  of  the  prominent  aneurisms 
seen  on  the  front  side  blast  wave. 

(v)  Due  to  the  high  temperature  in  the  forward  half  of  the  cavity  the 
flow  is  subsonic.  A  bow  shock  formed  about  a  small,  stationary  obstacle  to 
the  blast  wave  would  be  expected  to  dissipate  into  a  sound  wave  with  a  weaker 
density  jump  as  the  obstacle  is  engulfed  by  the  expanding  cavity. 
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igure  1.  Geometry  of  target,  laser  beam,  and  boundary  of  oblate  spheroi-dal 
coordinate  system  used  in  the  numerical  simulation.  The  system  is 
rotationally  symmetric  about  the  laser  beam  axis.  The  opening 
angle  for  the  blowoff  asymptotically  approaches  30°  far  from  the 
target.  The  front  and  back  boundary  are  each  1.1  cm  from  the 
target.  The  target  is  a  disc  of  radius  132  <j  and  thickness  J4.6  p. 
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3i.de  (o)  and  tack  side  ([])  blast  wave 
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The  position  of  the  front 
for  the  standard  model  as  a  function  of  time  measured  frcm 
peak  intensity  of  the  laser  pulse.  The  equal  mass  radii  for 
sides  are  al30  indicated.  The  solid  line  is  the  position  of 
front  side  blast  wave  frcm  the  experimental . 
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Figure  3.  The  energy  budget  for  the  standard  model  as  a  function  of  time 
measured  from  the  peak  intensity  of  the  laser  pulse.  /E^dt  is  the 
accumulated  energy  input  from  the  laser,  E^  is  the  instantaneous 
kinetic  energy  of  directed  motion,  E^  is  the  thermal  energy  of 
random  motion,  Ec  is  the  chemical  potential  energy,  and  /ER dt  is 
the  accumulated  radiated  energy. 
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gure  itd.  The  top  left  box  shows  the  rate  of  emitted  radiation  per  unit 
volume  (RAD  divided  by  1*10^)  and  the  rate  of  heating  per  unit 
volume  (HRAD)  due  to  absorption  of  the  radiation.  The  top  right 
box  breaks  down  the  emitted  radiation  into  recombination  'HR), 
line  (L),  and  bremstrahlung  (3).  The  bottom  left  box  shows  the 
electron  temperature  (TE)  and  an  effective  radiation  temperature 
(TRAD)  based  on  a  blackbody  flux.  The  bottom  right  box  shows  the 
smoothed  spectrum  emerging  from  the  back  side  (solid  line)  and 
front  3ide  (dotted  line).  The  simulation  time  of  6.66  nsec  is 
shortly  after  the  peak  intensity  of  the  laser  pulse  (5  nsec). 


Figure  5a.  Number  density  (NI  •  ions,  NE  ■  electrons,  NN  *  neutrals),  mass 
velocity,  pressure  (PI  -  ions,  PE  -  electrons,  Q  -  artificial 
viscosity)  and  temperature  (TI  -  ions,  TE  -  electrons)  for  model  B 


at  the  simulation  time  52.6  nsec. 


The  local  adiabatic  (Y  - 
position  for  four  different 
horizontal  dashed  line  is  at 


5/3)  Mach  number  as  a  function  of 
times  in  the  standard  model.  The 
the  value  unity. 
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1720  RANDOLPH  ROAD,  S.E. 

ALBUQUERQUE,  NM  87106 
01 CY  R.  STELLINGWEP.F 
01 CY  M.  ALME 
Old  L.  WRIGHT 

MITRE  CORP 

WESTGATE  RESEARCH  PARK 
1320  DOLLY  MADISON  BLVD 
MCLEAN,  VA  22101 

0 1 CY  ATTN  W.  HALL 
01 CY  ATTN  W.  FOSTER 

PACIFIC-SIERRA  RESEARCH  CORP 
1  a  3  0  SANTA  MONICA  BLVD. 

LOS  ANGELES,  CA  90025 

01 CY  ATTN  E.C.  FIELD,  JR. 

PENNSYLVANIA  STATE  UNIVERSITY 
IONOSPHERE  RESEARCH  LAB 
313  ELECTRICAL  ENGINEERING  EAST 
UNIVERSITY  PARK,  PA  16302 
(NO  CLASS  TO  THIS  ADDRESS) 

01 CY  ATTN  IONOSPHERIC  RESEARCH  LAB 

PHOTOMETRIES ,  INC. 

-  ARROW  DRIVE 
WOBURN,  MA  01801 

01 CY  ATTN  IRVING  L.  KOFSKY 


PHYSICAL  DYNAMICS,  INC. 

P.O.  BOX  10367 
OAKLAND.  CA  9*610 
ATTN  A.  THOMSON 

R  4  D  ASSOCIATES 
P.O.  BOX  9695 
MARINA  DEL  REY,  CA  90291 
01 CY  ATTN  FORREST  GILMORE 
01 CY  ATTN  WILLIAM  B.  WRIGHT,  JR 
0 1 CY  ATTN  WILLIAM  J.  KARZ AS 
0 1 C Y  ATTN  H.  ORY 
0 1 CY  ATTN  C.  MACDONALD 

RAND  CORPORATION,  THE 
15*50  COHASSET  STREET 
VAN  NUYS,  CA  91*06 

01 CY  ATTN  CULLEN  CRAIN 
0 1 C Y  ATTN  ED  5EDR0ZI AN 

RAYTHEON  CO. 

528  BOSTON  POST  ROAD 
SUDBURY,  MA  01776 

01 CY  ATTN  BARBARA  ADAMS 

RIVERSIDE  RESEARCH  INSTITUTE 
330  WEST  *  2nd  STREET 
NEW  YORK,  NY  10036 

0 1 C Y  ATTN  VINCE  TRAPANI 


SCIENCE  APPLICATIONS,  INC. 
1 1 50  PROSPECT  PLAZA 


JOLLA 

,  CA 

92037 

0 1  CY 

ATTN 

LEWIS  M.  LINSON 

01  CY 

ATTN 

DANIEL  A.  HAMLI 

01  CY 

ATTN 

E.  FRIEMAN 

01  CY 

ATTN 

E.A.  STRAKER 

Cl  CY 

ATTN 

CURTIS  A.  SMITH 

SCIENCE  APPLICATIONS,  INC 
1710  GOCDRIDGE  DR. 

MCLEAN,  VA  22102 
O’ CY  J.  COCKAYNE 
01 CY  E .  HYMAN 


PHYSICAL  DYNAMICS,  INC. 
P.O.  30X  3027 
3ELLEVUE,  WA  98009 

0 1 C Y  ATTN  E.J.  FREMCUW 
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SRI  INTERNATIONAL 
333  RAVENSWOOD  AVENUE 
MENLO  PARK ,  CA  9<*025 
0 1  C  Y  ATTN  J.  CASPER 
01  CY  ATTN  DONALD  NEILSON 
01 CY  ATTN  ALAN  BURNS 
01 CY  ATTN  G.  SMITH 
01 CY  ATTN  R.  TSUNODA 
01 CY  ATTN  DAVID  A.  JOHNSON 
01 CY  ATTN  HALTER  G.  CHESNUT 
01 CY  ATTN  CHARLES  L.  RINO 
01 CY  ATTN  HALTER  JAYE 
0 1 C Y  ATTN  J.  VICKREY 
01 CY  ATTN  RAY  L.  LEADA3RAND 
01 CY  ATTN  G.  CARPENTER 
01 CY  ATTN  G.  FRICE 
0 1 C Y  ATTN  R.  LIVINGSTON 
01 CY  ATTN  V.  GONZALES 
01 CY  ATTN  D.  MCDANIEL 

TECHNOLOGY  INTERNATIONAL  CORP 
75  HIGGINS  AVENUE 
BEDFORD,  MA  01730 

0 1 C Y  ATTN  H.P.  BOOUIST 

TRW  DEFENSE  A  SPACE  SYS  GROUP 
ONE  SPACE  PARK 
REDONDO  BEACH,  CA  90278 
0 1 C Y  ATTN  R.  K.  PLEBUCH 
01CY  ATTN  S.  ALTSCHULER 
01 CY  ATTN  D.  DEE 
0 1 C  Y  ATTN  0/  STOCKHELL 
SNTF/1575 

VISIDYNE 

SOUTH  3EDF0RD  STREET 
3u  RLINGTON ,  MA  0’803 
0 i C Y  ATTN  U.  REIDY 
C 1 C Y  ATTN  J.  CARPENTER 
0 1 C Y  ATTN  C.  HUMPHREY 

UNIVERSITY  OF  PITTSBURGH 
PITTSBURGH,  =A  *52’ 3 

01 CY  ATTN:  N.  ZA3USKY 

DIRECTOR  OF  RESEARCH 
l\  S.  NAVAL  ACADEMY 
ANNAPOLIS,  MD  21402 
02CY 

CODE  1220 
01  CY 
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